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One of the best motivated hypotheses in cosmology states that most of the matter in the universe is 
in the form of weakly-interacting massive particles that decoupled early in the history of the universe 
and cooled adiabatically to an extremely low temperature. Nevertheless, the hnite temperature and 
horizon scales at which these particles decoupled imprint generic signatures on their small scale den- 
sity fluctuations. We show that the previously recognized cut-off in the fluctuation power-spectrum 
due to free-streaming of particles at the thermal speed of decoupling, is supplemented by acoustic 
oscillations owing to the initial coupling between the cold dark matter (CDM) and the radiation 
field. The power-spectrum oscillations appear on the scale of the horizon at kinematic decoupling 
which corresponds to a mass scale of ~ 10~ 4 (Td/10 MeV) -3 Mq for a CDM decoupling temperature 
Td. The suppression of the power-spectrum on smaller scales by the acoustic oscillations is physi- 
cally independent from the free-streaming effect, although the two cut-off scales are coincidentally 
comparable for Td ~ 10 MeV and a particle mass of M ~ 100 GeV. The initial conditions for recent 
numerical simulations of the earliest and smallest objects to have formed in the universe, need to be 
modified accordingly. The smallest dark matter clumps may be detectable through 7-ray production 
from particle annihilation, through fluctuations in the event rate of direct detection experiments, 
or through their tidal gravitational effect on wide orbits of objects near the outer edge of the solar 
system. 

PACS numbers: 95.35.+d, 14.80.Ly, 98.35. Ce, 98.80.Cq 



Introduction. A broad range of observational data in- 
volving the dynamics of galaxies, the growth of large- 
scale structure, and the dynamics and nucleosynthesis of 
the universe as a whole, indicate the existence of dark 
matter with a mean cosmic mass density that is ~ 5 
times larger than the density of the baryonic matter 
[1,2]. The data is consistent with a dark matter com- 
posed of weakly-interacting, massive particles, that de- 
coupled early and adiabatically cooled to an extremely 
low temperature by the present time [1]. The Cold Dark 
Matter (CDM) has not been observed directly as of yet, 
although laboratory searches for particles from the dark 
halo of our own Milky- Way galaxy, have been able to 
restrict the allowed parameter space for these particles. 
Since an alternative more-radical interpretation of the 
dark matter phenomenology involves a modification of 
gravity [3], it is of prime importance to find direct fin- 
gerprints of the CDM particles. One such fingerprint 
involves the small-scale structure in the universe [4], on 
which we focus in this paper. 

Perhaps the most popular candidate for the CDM par- 
ticle is a Weakly Interacting Massive Particle (WIMP). 
The lightest supersymmetric particle (LSP) could be a 
WIMP (for a review see [1]). The CDM particle mass 
depends on free parameters in the particle physics model 
but typical values cover a range around M ~ 100 GeV 
(up to values close to a TeV). In many cases the LSP 
hypothesis will be tested at the Large Hadron Collider 
(e.g. [5]) or in direct detection experiments (e.g. [6]). 

The properties of the CDM particles affect their re- 



sponse to the small-scale primordial inhomogeneities pro- 
duced during cosmic inflation. The particle cross-section 
for scattering off standard model fermions sets the epoch 
of their thermal and kinematic decoupling from the cos- 
mic plasma (which is significantly later than the time 
when their abundance freezes-out at a temperature T ~ 
M). Thermal decoupling is defined as the time when the 
temperature of the CDM stops following that of the cos- 
mic plasma while kinematic decoupling is defined as the 
time when the bulk motion of the two species start to 
differ. For CDM the epochs of thermal and kinetic de- 
coupling coincide. They occur when the time it takes for 
collisions to change the momentum of the CDM particles 
equals the Hubble time. The particle mass determines 
the thermal spread in the speeds of CDM particles, which 
tends to smooth-out fluctuations on very small scales due 
to the free-streaming of particles after kinematic decou- 
pling [4,7]. Viscosity has a similar effect before the CDM 
fluid decouples from the cosmic radiation fluid [8] . An im- 
portant effect that has been previously ignored involves 
the memory the CDM fluid has of the acoustic oscilla- 
tions of the cosmic radiation fluid out of which it de- 
coupled. Here we consider the imprint of these acoustic 
oscillations on the small-scale power spectrum of density 
fluctuations in the universe. Other imprints of acoustic 
oscillations on much larger scales were identified recently 
in maps of the Cosmic Microwave Background (CMB) 
[2] , and the distribution of nearby galaxies [9] . 

Throughout this paper, we adopt the standard set of 
cosmological parameters [2] for a universe dominated by 
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cold dark matter and a cosmological constant (A-CDM). 

Formalism Kinematic decoupling of CDM occurs 
during the radiation-dominated era. For example, if 
the CDM is made of neutralinos with a particle mass 
of ~ 100 GeV, then kinematic decoupling occurs at a 
cosmic temperature of T<j <~ 10 MeV [8,10]. As long as 
Td <C 100 MeV, we may ignore the imprint of the QCD 
phase transition on the CDM power spectrum [11]. Over 
a short period of time during this transition, the sound 
speed of the plasma vanishes, resulting in a significant 
growth for perturbations with periods shorter than the 
length of time over which the sound speed vanishes. The 
transition occurs when the temperature of the cosmic 
plasma is <~ 100 — 200 MeV and lasts for a small fraction 
of the Hubble time. As a result, the induced modifica- 
tions are on scales smaller than those we are consider- 
ing here and the imprint of the QCD phase transition is 
washed-out by the effects we calculate. 

At early times the contribution of the dark matter to 
the energy density is negligible. Only at relatively late 
times when the cosmic temperature drops to values as 
low as ~ 1 eV, matter and radiation have comparable 
energy densities. As a result, the dynamics of the plasma 
at earlier times is virtually unaffected by the presence of 
the dark matter particles. In this limit, the dynamics of 
the radiation determines the gravitational potential and 
the dark matter just responds to that potential. We will 
use this simplification to obtain analytic estimates for the 
behavior of the dark matter transfer function. 

The primordial inflationary fluctuations lead to acous- 
tic modes in the radiation fluid during this era. The 
interaction rate of the particles in the plasma is so high 
that we can consider the plasma as a perfect fluid down 
to a comoving scale, 
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where % = f* d dt/a(t) is the conformal time (i.e. the 
comoving size of the horizon) at the time of CDM decou- 
pling, td; o is the scattering cross section and n is the rel- 
evant particle density. (Throughout the paper we set the 
speed of light and Planck's constant to unity.) The damp- 
ing scale depends on the species being considered and its 
contribution to the energy density and is the largest for 
neutrinos which are only coupled through weak interac- 
tions. In that case N ~ {T/T»f where T% ~ 1 MeV is 
the temperature of neutrino decoupling. At the time of 
CDM decoupling N ~ M/T d ~ 10 4 for the rest of the 
plasma, where M is the mass of the CDM particle. In 
this paper we will consider modes of wavelength larger 
than A/, and so we neglect the effect of radiation diffu- 
sion damping and treat the plasma (without the CDM) 
as a perfect fluid. 

The equations of motion for a perfect fluid during the 
radiation era can be solved analytically. We will use that 
solution here, following the notation of Rcf. [12]. As usual 



we Fourier decompose fluctuations and study the behav- 
ior of each Fourier component separately. For a mode of 
comoving wavenumbcr k in Newtonian gauge, the gravi- 
tational potential fluctuations are given by: 



$ = 3$ E 



sin^Lorj) — lot] cos(lot]) 



(2) 



where w = k/^/3 is the frequency of a mode and $ p is its 
primordial amplitude in the limit rj — > 0. In this paper 
we use conformal time r] = J dt/a(t) with a(t) oc t 1 / 2 
during the radiation-dominated era. The monopole ©o 
and dipole ©i of the photon distribution can be written 
in terms of the gravitational potential as: 
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where x = kr\ and a prime denotes a derivative with 
respect to x. 

The solutions in equations (2) and (3) assume that 
both the sound speed and the number of relativistic de- 
grees of freedom are constant over time. As a result of 
the QCD phase transition and of various particles be- 
coming non-rclativistic, both of these assumptions are 
not strictly correct. The vanishing sound speed during 
the QCD phase transition provides the most dramatic ef- 
fect, but its imprint is on scales smaller than the ones we 
consider here because the transition occurs at a signifi- 
cantly higher temperature and only lasts for a fraction of 
the Hubble time [11]. 

Before the dark matter decouples kinematically we will 
treat it as a fluid which can exchange momentum with 
the plasma through particle collisions. At early times, 
the CDM fluid follows the motion of the plasma and is 
involved in its acoustic oscillations. The continuity and 
momentum equations for the CDM can be written as: 

5 C + 6 C = 3$ 

9 C + U c = k 2 c 2 s 5 c - k 2 a c - fc 2 $ + r-^Gi - 6 C ) (4) 

where a dot denotes an ry-derivative, 8 C is the dark mat- 
ter density perturbation, 8 C is the divergence of the 
dark matter velocity field and a c denotes the anisotropic 
stress. In writing these equations we have followed Ref. 
[13]. The term t c _1 (@i — 9 C ) encodes the transfer of mo- 
mentun between the radiation and CDM fluids and t^ 1 
provides the collisional rate of momentum transfer, 
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with n being the number density of particles with which 
the dark matter is interacting, a{T) the average cross sec- 
tion for interaction and M the mass of the dark matter 
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particle. The relevant scattering partners are the stan- 
dard model leptons which have thermal abundances. For 
detailed expressions of the cross section in the case of su- 
persymmetric (SUSY) dark matter, see Refs. [10,7]. For 
our purpose, it is sufficient to specify the typical size of 
the cross section and its scaling with cosmic time, 
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where the coupling mass M a is of the order of the weak- 
interaction scale (~ 100 GeV) for SUSY dark matter. 
This equation should be taken as the definition of M a , 
as it encodes all the uncertainties in the details of the 
particle physics model into a single parameter. The tem- 
perature dependance of the averaged cross section is a 
result of the available phase space. Our results are quite 
insensitive to the details other than through the decou- 
pling time. Equating r ( T 1 /a to the Hubble expansion rate 
gives the temperature of kinematic decoupling: 
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The term k 2 c 2 S c in Eq. (4) results from the pressure 
gradient force and c s is the dark matter sound speed. In 
the tight coupling limit, r c <C H^ 1 we find that c 2 s w 
f c T/M and that the shear term is k 2 a c w f v c 2 T c 9 c . Here 
/„ and f c are constant factors of order unity. We will 
find that both these terms make a small difference on the 
scales of interest, so their precise value is unimportant. 

By combining both equations in (4) into a single equa- 
tion for S c we get 
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where Xd = krj d and r\ d denotes the time of kinematic 
decoupling which can be expressed in terms of the de- 
coupling temperature as, 



rj d = 2t d (l + z d ) « « 10 pc / T( ' 



T T d ~ r ~ \10 MeV, 

cx M^M" 1 / 4 , (9) 

with To = 2.7K being the present-day CMB temperature 
and z d being the redshift at kinematic decoupling. We 
have also introduced the source function, 
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For iCi,;, the dark matter sound speed is given by 

x d 



c s (x) = c s (x d )'- 
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where c 2 (x d ) is the dark matter sound speed at kinematic 
decoupling (in units of the speed of light), 
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In writing (11) we have assumed that prior to decoupling 
the temperature of the dark matter follows that of the 
plasma. For the viscosity term we have, 



F v (x)=f v c 2 s (x d )x 2 (^y 



(13) 



Free streaming after kinematic decoupling In the limit 
of the collision rate being much slower than the Hub- 
ble expansion, the CDM is decoupled and the evolution 
of its perturbations is obtained by solving a Boltzman 
equation: 



df | dxj df | dq l df = Q 
drj drj dxi dr\ dqi 



(14) 



where f(x,q,r)) is the distribution function which de- 
pends on time, position and comoving momentum q. The 
comoving momentum 3-components are dxi/drj = qi/a. 
We use the Boltzman equation to find the evolution of 
modes that are well inside the horizon with x » 1. In 
the radiation era, the gravitational potential decays after 
horizon crossing (see Eq. 2). In this limit the comoving 
momentum remains constant, dqi /drj = and the Boltz- 
man equation becomes, 



^L + ^iK = 

dt] a dxi 

We consider a single Fourier mode and write / as, 

f(x,q,r ] ) = f Q (q)[l + 8 F (q,r ] )e^% 
where fo(q) is the unperturbed distribution, 
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where iicdm and Tcdm are the present-day density and 
temperature of the dark matter. 

Our approach is to solve the Boltzman equation with 
initial conditions given by the fluid solution at a time 
?7» (which will depend on fc). The simplified Boltzman 
equation can be easily solved to give 5f(<1, tj) as a function 
of the initial conditions 5f{(I,ii*), 
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The CDM overdensity S c can then be expressed in 
terms of the perturbation in the distribution function as, 



3 



"CDM 



d 3 q fo(q) S F (q,r])- 



(19) 



We can use (18) to obtain the evolution of 5 C in terms of 
its value at 77*, 
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where fc7 2 = ^/ (Td/M)rjd. The exponential term is re- 
sponsible for the damping of perturbations as a result 
of free streaming and the dispersion of the CDM parti- 
cles after they decouple from the plasma. The above ex- 
pression is only valid during the radiation era. The free 
streaming scale is simply given by J dt(v/a) cx J dta~ 2 
which grows logarithmically during the radiation era as 
in equation (20) but stops growing in the matter era when 
a oc t 2 / 3 . 

Equation (20) can be used to show that even during the 
free streaming epoch, S c satisfies equation (8) but with a 
modified sound speed and viscous term. For x 3> Xd one 
should use, 
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(21) 



The differences between the above scalings and those dur- 
ing the tight coupling regime are a result of the fact that 
the dark matter temperature stops following the plasma 
temperature but rather scales as a~ 2 after thermal decou- 
pling, which coincides with the kinematic decoupling. We 
ignore the effects of heat transfer during the fluid stage 
of the CDM because its temperature is controlled by the 
much larger heat reservoir of the radiation-dominated 
plasma at that stage. 

To obtain the transfer function we solve the dark mat- 
ter fluid equation until decoupling and then evolve the 
overdensity using equation (20) up to the time of matter- 
radiation equality. In practice, we use the fluid equations 
up to x* = 10 max.(xd, 10) so as to switch into the free 
streaming solution well after the gravitational potential 
has decayed. In the fluid equations, we smoothly match 
the sound speed and viscosity terms at x = x d . As men- 
tioned earlier, because c s (xd) is so small and we are in- 
terested in modes that are comparable to the size of the 
horizon at decoupling, i.e. Xd ~ few, both the dark mat- 
ter sound speed and the associated viscosity play only a 
minor role, and our simplified treatment is adequate. 

In Figure 1 we illustrate the time evolution of modes 
during decoupling for a variety of k values. The situation 
is clear. Modes that enter the horizon before kinematic 
decoupling oscillate with the radiation fluid. This be- 
havior has two important effects. In the absence of the 




FIG. 1. The normalized amplitude of CDM fluctuations 
<5/$p for a variety of modes with comoving wavenumbers 
log(krjd) — (0, 1/3, 2/3, 1, 4/3, 5/3, 2) as a function of x = krj, 
where 77 = J dt/a(t) is the conformal time coordinate. The 
dashed line shows the temperature monopole 3#o and the dot- 
ted curve shows the evolution of a mode that is uncoupled to 
the cosmic plasma. 



coupling, modes receive a "kick" by the source term S(x) 
as they cross the horizon. After that they grow logarith- 
mically. In our case, modes that entered the horizon be- 
fore kinematic decoupling follow the plasma oscillations 
and thus miss out on both the horizon "kick" and the 
beginning of the logarithmic growth. Second, the decou- 
pling from the radiation fluid is not instantaneous and 
this acts to further damp the amplitude of modes with 
id > 1. This effect can be understood as follows. Once 
the oscillation frequency of the mode becomes high com- 
pared to the scattering rate, the coupling to the plasma 
effectively damps the mode. In that limit one can replace 
the forcing term 9g by its average value, which is close 
to zero. Thus in this regime, the scattering is forcing the 
amplitude of the dark matter oscillations to zero. After 
kinematic decoupling the modes again grow logarithmi- 
cally but from a very reduced amplitude. The coupling 
with the plasma induces both oscillations and damping of 
modes that entered the horizon before kinematic decou- 
pling. This damping is different from the free streaming 
damping that occurs after kinematic decoupling. 

In Figure 2 we show the resulting transfer function of 
the CDM overdensity. The transfer function is defined 
as the ratio between the CDM density perturbation am- 
plitude S c when the effect of the coupling to the plasma 
is included and the same quantity in a model where the 
CDM is a perfect fluid down to arbitrarily small scales 
(thus, the power spectrum is obtained by multiplying the 
standard result by the square of the transfer function). 
This function shows both the oscillations and the damp- 
ing signature mentioned above. The peaks occur at mul- 
tipoles of the horizon scale at decoupling, 
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FIG. 2. Transfer function of the CDM density pertur- 
bation amplitude (normalized by the primordial amplitude 
from inflation). We show two cases: (i) T d /M = 1CT 4 and 
T d /T cq = 10 7 ; (ii) T d /M = 1(T 5 and T d /T oq = 10 7 . In each 
case the oscillatory curve is our result and the other curve is 
the free-streaming only result that was derived previously in 
the literature [4,7,8]. 
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This same scale determines the "oscillation" damping. 
The free streaming damping scale is, 



ridCd(r)d) \n(r) eq /i] d ) oc 
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(23) 



where T cq is the temperature at matter radiation equal- 
ity, T eq w 1 eV. The free streaming scale is parametri- 
cally different from the "oscillation" damping scale. How- 
ever for our fiducial choice of parameters for the CDM 
particle they roughly coincide. 

The vanishing of the sound speed during the QCD 
phase transition amplifies perturbations which have 
luAtjqcd > 1, where Atjqcd is the duration of the tran- 
sition. In Figure 2 the affected modes are those with 
Xd = kr] d > a/3 (tiqcd/Avqcd) (Vd/VQCD)- Typical 
values of rj d /r) QCD — 10-15 and t\qcdI Atjqcd ~ 3 - 10 
relate this condition to modes with x d > 50 — 260. Thus 
the affected scales are severely damped by the effects con- 
sidered in this paper. 

Finally we want to stress the fact that the damping 
scale is significantly smaller than the scales observed di- 
rectly in the Cosmic Microwave Background or through 
large scale structure surveys. For example, the ratio of 
the damping scale to the scale that entered the horizon 
at equality is rj d /rj eq ~ T eq /T d ~ 10~ 7 and to our present 
horizon %/% ~ (T eg T ) 1/2 /T d ~ 10~ 9 . In the context 
of inflation, these scales were created 16 and 20 e-folds 
apart. Given the large extrapolation, one could certainly 
imagine that a change in the spectrum could alter the 
shape of the power spectrum around the damping scale. 



However, for smooth inflaton potentials with small de- 
partures from scale invariance this is not likely to be the 
case. On scales much smaller than the horizon at matter 
radiation equality, the spectrum of perturbations density 
before the effects of the damping are included is approx- 
imately, 



A 2 (fc) oc exp 



(n - 1) ln(krj eq ) + \<? \n(kr] eq ) 2 + 



x ln 2 (fc?7 e9 /8) 



(24) 



where the first term encodes the shape of the primordial 
spectrum and the second the transfer function. Primor- 
dial departures from scale invariance are encoded in the 
slope n and its running a. The effective slope at scale k 
is then, 



d In A 2 
d In k 



= (n - 1) + ala(kr}eq) + 



\n(k Veq /8) ' 



(25) 



For typical values of (n — 1) ~ 1/60 and a ~ 1/60 2 
the slope is still positive at k ~ rQ 1 , so the cut-off in the 
power will come from the effects we calculate rather than 
from the shape of the primordial spectrum. However 
given the large extrapolation in scale, one should keep in 
mind the possibility of significant effects resulting from 
the mechanisms that generates the density perturbations. 

Implications We have found that acoustic oscilla- 
tions, a relic from the epoch when the dark matter cou- 
pled to the cosmic radiation fluid, truncate the CDM 
power spectrum on a comoving scale larger than effects 
considered before, such as free-streaming and viscosity 
[4,7,8]. For SUSY dark matter, the minimum mass of 
dark matter clumps that form in the universe is there- 
fore increased by more than an order of magnitude to a 
value of * 

47r ( t 

A/cut = — -, "M/Ocrit 

3 \k cnt J 



10- 



10 MeV J 



M , 



(26) 



where p crit = (Hq/8ttG) = 9 x 10 -30 g cm -3 is the crit- 
ical density today, and SIm is the matter density for the 
concordance cosmological model [2], We define the cut- 
off wavenumber fc cu t as the point where the transfer func- 
tion first drops to a fraction 1/e of its value at k — > 0. 
This corresponds to fc cut rs 3.3 rj^ 1 . 

Recent numerical simulations [15,16] of the earliest and 
smallest objects to have formed in the universe [17], need 



*Our definition of the cut-off mass follows the convention of 
the Jeans mass, which is defined as the mass enclosed within 
a sphere of radius Aj/2 where Aj = 2n/kj is the Jeans wave- 
length [14]. 



5 



to be redone for the modified power spectrum that we cal- 
culated in this paper. Although it is difficult to forecast 
the effects of the acoustic oscillations through the stan- 
dard Press-Schechter formalism [18], it likely that the 
results of such simulations will be qualitatively the same 
as before except that the smallest clumps would have a 
mass larger than before (as given by Eq. 26). 

Potentially, there are several observational signatures 
of the smallest CDM clumps. As pointed out in the liter- 
ature [15,19], the smallest CDM clumps could produce 7- 
rays through dark-matter annihilation in their inner den- 
sity cusps, with a flux in excess of that from nearby dwarf 
galaxies. If a substantial fraction of the Milky Way halo 
is composed of CDM clumps with a mass ~ 10 _4 M©, the 
nearest clump is expected to be at a distance of <~ 4x 10 17 
cm. Given that the characteristic speed of such clumps 
is a few hundred km s _1 , the 7-ray flux would therefore 
show temporal variations on the relatively long timescale 
of a thousand years. Passage of clumps through the solar 
system should also induce fluctuations in the detection 
rate of CDM particles in direct search experiments. 

Other observational effects have rather limited 
prospects for detectability. Because of their relatively 
low-mass and large size (~ 10 17 cm), the CDM clumps 
are too diffuse to produce any gravitational lensing sig- 
natures (including /emto-lensing [20] ) , even at cosmolog- 
ical distances. They could, however, have an effect on 
the time delay of the signal from transient cosmological 
sources. We illustrate this effect by considering a charac- 
teristic CDM clump with a size ~ 10 17 cm and a mass of 
10~ 4 M© (corresponding to formation rcdshift ~ 50). If 
a major fraction of the intergalactic CDM is in the form 
of these clumps (amounting to a mean mass density of 
~3x 10~ 30 g cm -3 ), then the mean separation between 
clumps is ~ 6 x 10 19 cm= 20pc. The gravitational poten- 
tial would then fluctuate on this length scale everywhere 
in the universe with an amplitude of ($/c ) ~2x 10~ 19 . 
Now, consider a compact source at cosmological distances 
with a size comparable to the distance between clumps, 
~ lOpc. Photons emerging from different regions of the 
source will encounter different gravitational potential de- 
lays (the so-called Shapiro time delay [21]) along their 
path due to the fact that they traverse CDM clumps at 
different impact parameters. For a source that is much 
larger than ~ 20pc, the effect of multiple CDM clumps 
would average out because the source would sample the 
full distribution of impact parameters with small statis- 
tical fluctuations around the mean time delay. The mean 
Shapiro delay from the network of CDM clumps along the 
line-of-sight to a cosmological source is ~ ($/c 2 )t# <~ 0.1 
sec, where tn is the light travel time across the observable 
universe. The total number of clumps along the line-of- 
sight is N ~ (4Gpc/20pc) - 2 x 10 8 . This implies that 
the variance in time delays across a source of size ~ lOpc 
would be suppressed by a factor of ~ l/y/N <~ 10~ 4 rel- 
ative to the mean time delay, i.e. it would have a mag- 



nitude ~ 10 /isec. 7-ray bursts may possess variability 
on such a timescale (which would have been smoothed 
out if their source had a size ~ lOpc), but their esti- 
mated source size is R < O.lpc [22] and so the tempo- 
ral smoothing effect for them is smaller by another fac- 
tor > (20/0.1) = 2 x 10 2 , which leads to a variance in 
the delay timescale of less than a microsecond. Unfortu- 
nately this time delay appears unmeasurable, since the 
variability of any source of apparent dimension R would 
be expected to be smoothed-out on a timescale, ~ R/jc, 
where 7 is the Lorentz factor of the emitting material 
(which is < 10 3 for 7-ray bursts). 

The smallest CDM clumps should not affect the inter- 
galactic baryons which have a much larger Jeans mass 
[17]. However, once objects above ~ 10 6 M Q start to col- 
lapse at redshifts z < 30, the baryons would be able to 
cool inside of them via molecular hydrogen transitions 
and the interior baryonic Jeans mass would drop. The 
existence of dark matter clumps could then seed the for- 
mation of the first stars inside these objects [23]. 

Finally, we note that the smallest CDM clumps may 
have a dynamical effect on wide orbits of test particles 
near the outer edge the solar system after the effects of 
the known planets have been modeled to a high precision 
[24]. 
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